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ABSTRACT 

We present relations of the black hole mass and the optical luminosity with the velocity dispersion 
and the luminosity of the [Ne V] and the [O IV] high-ionization lines in the mid-infrared (MIR) for 28 
reverberation-mapped active galactic nuclei. We used high-resolution Spitzer Infrared Spectrograph 
and Infrared Space Observatory Short Wavelength Spectrometer data to fit the profiles of these MIR 
emission lines that originate from the narrow- line region of the nucleus. We find that the lines are 
often resolved and that the velocity dispersion of [Ne V] and [O IV] follows a relation similar to that 
between the black hole mass and the bulge stellar velocity dispersion found for local galaxies. The 
luminosity of the [Ne V] and the [O IV] lines in these sources is correlated with that of the optical 
5100 A continuum and with the black hole mass. Our results provide a means to derive black hole 
properties in various types of active galactic nuclei, including highly obscured systems. 
Subject headings: infrared: galaxies — galaxies: active — galaxies: nuclei — galaxies: Seyferts — 
galaxies: kinematics and dynamics — quasars: emission lines 



1. INTRODUCTION 

Following the relation between the black hole (BH) 
mass, Mbh, and the bulge stellar velocity dispersion, 
cr,, (Ferrarese & Merritt 2000; Gebhardt et al. 2000; 
Tremaine et al. 2002), a similar relation was found for 
the velocity dispersion a of the [O III] 5007 A line (Nel- 
son 2000; Greene & Ho 2005) that originates from the 
narrow-line region (NLR) gas of the active galactic nu- 
cleus (AGN). The NLR gas is thought to be mostly grav- 
itationally bound to the bulge in high-luminosity AGNs 
(Nelson & Whittle 1996; Greene & Ho 2005; Laor 2007), 
unlike the broad-line region gas that is virialized due its 
proximity to the BH (Peterson & Wandel 1999). This 
relation is useful for systems in which the stellar absorp- 
tion lines cannot be observed, e.g., because of dilution 
of the stellar light by the AGN continuum. In addition 
to this relation, Mbh and the luminosity L of the 5100 
A optical continuum, Lopt, are correlated in AGNs with 
BH measurements (Kaspi et al. 2000). 

There are two main advantages to expanding such re- 
lations in the mid-infrared (MIR). The first is that the 
Ne V and O IV ions emitting at 14.32 and 25.89 /im can- 
not be easily excited by star-forming regions since they 
have ionization potentials x of 97.12 and 54.93 eV. The 
second reason is the low obscuration, which allows for 
the results to be apphed to type 2 AGNs. In this Letter, 
we investigate for relations between Mbh and Lopt with 
MIR NLR line velocity dispersions and luminosities. 

2. THE SAMPLE AND THE DATA REDUCTION 
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Our sources were selected from a sample of 35 local 
AGNs that have Mbh measurements from reverbera- 
tion experiments (Peterson et al. 2004). Archival high- 
resolution Spitzer Infrared Spectrograph (IRS) data are 
available for 28 reverberation-mapped AGNs. Of these, 
eight were previously observed with the Infrared Space 
Observatory Short Wavelength Spectrometer (SWS) in 
high resolution (Sturm et al. 2002). 

To reduce the IRS data we first averaged all images 
at the same nod position, rejecting pixels that deviated 
from the median image by more than 5 times its stan- 
dard deviation. We computed the uncertainty image and 
then created the bad-pixel mask and extracted the spec- 
trum of each nod position using the Spitzer data reduc- 
tion packages IRSCLEAN and SPICE. For the spectral 
extraction, we assumed a point-source spatial distribu- 
tion for the high-ionization lines, corresponding to an 
unresolved NLR. When present, we removed the weak 
[CI II] and [Fe II] lines at 14.37 and 25.99 ^m, 1000 and 
1200 kms-i longward of [Ne V] and [OIV]. We fit- 
ted the [Ne V] and [O IV] lines and their underlying 
continua using Gaussian and second-order polynomial 
functions, respectively, which was a good approximation 
for most cases. We considered resolved all lines with 
f-m — Cm, > cTi + Ci, whcrc (Ti is the instrumental resolu- 
tion i? at a given wavelength divided by 2.35, is the 
error of Ui, and am is the measured velocity dispersion. 
The error of am, f-m, is computed as (e^^^ -I- e^)°-^. The 
quantity tstd is the standard deviation of the am val- 
ues that are obtained when using different polynomial 
coefficients to describe the local continuum. It encap- 
sulates uncertainties in the continuum level and slope. 
The instrumental resolution and its error were measured 
as a function of wavelength using standard IRS calibra- 
tion targets (P Cygni, HD 190429, HD 174638). The 
average value of R in the 14.0—18.0 fMn range is 494±59 
kms~^. Between 25.0 and 34.2 /itm, it is equal to 503±63 
kms~^. The intrinsic velocity dispersions of the resolved 
lines were computed as (cr^ ~ ^fY''^ and are presented in 
Table 1 . Examples of lines resolved with IRS are given in 
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Figure 1. To complement the Spitzer data and to com- 
pute the widths of unresolved lines, we used SWS spectra 
presented in Sturm et al. (2002). SWS has a higher reso- 
lution {R ~160 kms-i at 14.32 fim and R ~230 kms-^ 
at 25.89 /xm) but lower sensitivity than IRS. It detected 
only the brightest [Ne V] and [O IV] lines, the velocity 
dispersions of which are presented in Table 1. To fit the 
relations between the MIR line and the BH properties, 
we used the FITEXY algorithm that takes into account 
the error bars in both axes (Press et al. 1992) and that 
treats all variables symmetrically (Tremaine et al. 2002). 

3. RESULTS 

3.1. Correlations between Mbh and the velocity 

dispersion of [Ne V] and [O IV]. 

The Mbh — c relation of the MIR high- ionization lines 
is presented in Figure 2. Its best-fit solution is 

(^) = «1 + X log (2001^) ' (1) 

with parameter values corresponding to ai = 7.86 ± 0.09 
and Pi = 4.31 ± 0.63 for [Ne V] and = 8.04 ± 0.08 
and I3i = 4.89 ± 0.63 for [O IV]. The scatter values of 
the relation along the Mbh axis are 0.55 and 0.51 dex 
for [Ne V] and [O IV] , respectively. Its Spearman rank 
coefficients are 0.66 and 0.77, assigning probabilities p 
of 5.0x10^'^ and 8.3x10^'' to the null hypothesis that 
the variables are correlated by chance. We excluded all 
sources with Mbh< 10^ Mq from this fit since the IRS 
resolution is insufficient to provide a measurements on 
the left-hand side of the relation and since no SWS mea- 
surements are available for such AfBH values. Moreover, 
IC4329A is not used in the fit since its light curves are 
not indicative of the presence of a BH (Peterson et al. 
2004). Given that the best-fit slope agrees within the 
errors with that of the Tremaine et al. (2002) formula, 
the latter can be alternatively used to derive Mbh from 
(T[Ng VI ^iid (T[o IV] ■ If we fix the slope to /3i — 4.02 and 
solve for qi using the IDL routine MPFIT, we find that 
ai = 7.90 ± 0.25 and ai = 8.11 ± 0.26 for [Ne V] and 
[O IV]. 

3.2. Correlations between Lopt; Mbh> and the 

luminosity of [Ne V] and [O IV]. 

Following the suggestion that the luminosity of [O III] 
can be used as a proxy of Lopt (Heckman et al. 2004) , we 
find that L[Ne v] and L[o iv] are related to the 5100 A 
continuum luminosity (Fig. 3) as 

log ( w^tr^L-^ ) = + /?2 X log ( ) , (2) 

and they are computed for Hq=70 kms^^ Mpc^^, 
^^m=0.3, and ritotai=l- The best-fit parameter values 
are a2 = -0.04 ± 0.01 and /?2 = 1.17 ± 0.02 for [Ne V] 
and a2 = -0.55 ± 0.01 and [^2 = 1-15 ± 0.01 for [O IV]. 
The scatter values of the relation are 0.46 and 0.47 dex 
and its Spearman rank coefficients are 0.85 and 0.81 for 
[Ne V] and [O IV] , respectively, corresponding to proba- 
bilities p of 3.7x10"^ and 7.5x10"^. The relations con- 
tinue to hold when plotting the 5100 A continuum versus 
the MIR line fiuxes (Fig. 3). In the flux-flux diagram, 
the only outlier is PG 1226+023 (3C 273), which is the 
brightest Palomar-Green QSO (Schmidt & Green 1983). 
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Fig. 1. Examples of lines resolved with IRS. The [Ne V] pro- 
file is given for PG 2130+099 (solid line, diamonds) and Mrk279 
(dashed line, stars). The resolution element is plotted as a dotted 
line. 



A correlation similar to that found in the optical be- 
tween Mbh and Lopt exists in the MIR between Mbh 
and L[Ne v] or ^[o iv] (Fig- 4). Its best fit is 

l°g (^) = «3 + /33 X log ( >L%'L- ) > (3) 

with aa = 8.09 ± 0.02 and /Jg = 0.74 ± 0.03 for [Ne V], 
and as = 7.72 ± 0.02 and /^s = 0.77 ± 0.03 for [O IV]. 

Its scatter values are 0.46 and 0.50 dex and its Spearman 
rank coefficients are 0.76 and 0.69 with p = 2.5 x 10~^ 
and p = 1.3 X lO^'' for [Ne V] and [O IV], respectively. 
We note that upper limits have not been used for the fit. 

4. DISCUSSION: ORIGIN AND CONDITIONS FOR THE 
APPLICATION OF THE RELATIONS 

The relations between CT[no v] ^ (^[o iv] and Mbh indi- 
cate that the NLR gas kinematics are primarily deter- 
mined by the potential of the bulge, as it is also believed 
based on the [O III] line profiles (Whittle 1992; Nelson 
& Whittle 1996; Greene & Ho 2005). However, a^-^e v] 
and cr[o iv] are on average 67 and 62 kms"-"- higher than 
C7*, and 51 and 32 kms~^ higher than cr[o iii]- Such dis- 
crepancies could be attributed to an increase of the line 
width with increasing which is often found for optical 
NLR fines (e.g. Osterbrock 1981; Whittle 1985b). It 
is possible that ions with high ionization potentials have 
high velocity dispersions because they are located in NLR 
clouds that are close to the BH sphere of influence. Since 
the O IV and, mostly, the Ne V ions are predominantly 
excited by AGNs, the MIR Mbh — c relation can be ap- 
plied even in systems that undergo starbursts. It can 
also be applied in type 2 AGNs since the lines suffer lit- 
tle from extinction. Moreover, cr[Ne v] and (T[q iv] can 
be used as surrogates for cr* in environments where mea- 
suring (7* is hard, such as in bright QSOs. However, the 
relation does not necessarily hold for AGNs with strong 
winds and jets such as luminous radio sources (Whittle 
1992), and for AGNs with high Eddington rates, ryEdd 
(Greene & Ho 2005). Such systems can be recognized 
by the asvmmctric wings or the high kurtosis of their line 
profiles (Whittle 1985a; 1992). 

The luminosities of the broad lines (Baldwin 1977) 
and the narrow lines (Fig. 3) scale with iopt and with 
the bolometric AGN luminosity, Lboi, since they depend 
on the absolute accretion rate onto the central object. 
Specifically, iboi is equal to XoptLopt, a;[Ne v] -^[Ne v]. 



High-ionization MIR lines as AGN black hole mass and bolometric luminosity indicators 
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Fig. 2. — The relation between Mbh and the velocity dispersion of the MIR NLR lines [Ne V] (left, circles) and [O IV] (right, triangles). 
The solid line corresponds to the best-fit solution and the open symbols correspond to sources that have been excluded from the fit. For 
comparison, we overplot the Tremaine et al. (2002) relation (dashed line) and the stellar. Call, velocity dispersion measurements of the 
reverberation-mapped AGNs (stars; Onken et al. 2004; Nelson et al. 2004). The crosses correspond to cttq ttti. 
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Fig. 3. — The luminosity (left) and fiux (right) of the optical continuum scale with those of the MIR NLR lines. The best-fit solution to 
the [Ne V] data points (circles) is plotted with a dashed line and that to the [O IV] data points (triangles) is plotted with a dash-dotted 
line. The open symbols correspond to limits, not used for the fit. 
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Fig. 4. — The relation between Mbh and the luminosity of the 
MIR NLR lines. The symbols are as in Figure 3. The dotted lines 
correspond to a fixed Eddington rate of 0.01, 0.1, and 1 from top 
to bottom. They were horizontally shifted by a factor of 1000 to 
overlap with the luminosity range of the MIR NLR lines. 



and a;[o iv] -^[o iv]: where a; is a wavelength-dependent 
bolometric correction factor. For a;opt=9 (Kaspi et al. 
2000; ]VIarconi et al. 2004) and for the median values 
of Lo-pt/ i[No V] and -Lopt/ L^Q iv] in our sample, we find 
that xjNe V] =13000 and x^q iv]=4500. That the relation 
does not have a slope of unity could indicate that X[nc v] 
and a;[o iv] depend on L in a manner different than .Topt 
does. Netzer et al (2006) found that the equivalent width 
of [O III] decreases with increasing Lopt- This behaviour 



could be attributed to a different geometric distribution 
of NLR clouds around the AGN at different luminosities. 

Whether the [Ne V] or [O IV] luminosity of a source 
can be used to derive its Mbh partly depends on its Ed- 
dington rate. To illustrate how changes of T^Edd affect the 
relation between the luminosity of a MIR line and AfeH, 
we overplot lines of constant T/Edd in Figure 4. The more 
quiescent a source is, the more its position is shifted to 
the top left corner of this diagram. Equation (3) is valid 
within the ?7Edd range 0.003—0.6 that our sources span, 
and implies that the mass of a BH determines its lumi- 
nosity output. Its advantage is that it can be applied 
to sources with obscured optical continua or underesti- 
mated [O III] luminosities, such as type 2 AGNs (Netzer 
et al 2006). However its scatter is likely to be larger 
than that presented in § 3.2 since the Eddington rates 
of the reverberation-mapped AGNs are not necessarily 
representative of those of all local AGNs. 

We conclude that the MIR NLR gas kinematics trace 
Mbh in local reverberation-mapped AGNs in a manner 
similar to stellar kinematics. The [Ne V] and [O IV] line 
widths can be used to estimate Mbh- The calibration of 
-^[Nc V] and L[o iv] to Lopt provides a new method to 
compute bolometric luminosities, black hole masses, and 
Eddington rates in various types of AGNs, including 
highly obscured systems, albeit with a large scatter. 

This work was based on observations made with 
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TABLE 1 

Black hole and MIR line properties for reverberation-mapped AGNs with high-resolution MIR spectra 



Galaxy 


z 


Mbh 






<^[0 IIIJ 




''"[Ne V] 


P\o IV\ 


^■[O IV] 






(lO'' Mq) 


(10*4 ergs s-i) 


(km s-l) 


(km s-^) 


(10-21 W cm-2) 


(km s-l^) 


(10-21 W cm-2) 


(km s-i) 


3C 120 


0.03301 


5.55 ± 2.70 


1.479 ± 0.272 


162 ± 24 




18.06 ± 0.23 


111 ± le'' 


120.2 ± 0.8 


114 ± 9=- 


Fairall 9 


0.04702 


25.5 ± 5.6 


1.778 ± 0.205 




181 


2.60 ± 0.33 


208 ± 36 


6.01 ± 0.28 


215 ± 39 


IC 4329A 


0.01605 


0.99 ± 1.49 


0.208 ± 0.024 


122 ± 13 


234 


23.63 ± 1.22 


304 ± 44 


100.6 + 1.1 


278 ± 30 


Mrk 279 


0.03045 


3.49 ± 0.92 


0.756 ± 0.087 


197 ± 12 


247 


3.28 ± 0.23 


247 ± 35 


10.82 ± 0.29 




Mrk 509 


0.0344 


14.3 ± 1.2 


1.905 ± 0.351 




221 


6.29 ± 0.26 




28.73 ± 0.29 


133 ± 34=' 


Mrk 590 


0.02638 


4.75 ± 0.74 


0.646 ± 0.134 


192 ± 10 


170 


1.01 ± 0.16 




3.39 ± 0.17 


197 ± 30 


NGC 3227 


0.00386 


4.22 ± 2.14 


0.024 ± 0.002 


134 ± 6 


206 


23.13 ± 0.35 


212 ± 34 


65.37 ± 0.59 


210 ± 30 


NGC 3783 


0.00973 


2.98 ± 0.54 


0.182 ± 0.017 


95 ± 10 


98 


15.24 ± 1.00 


133 ± 21* 


39.78 ± 0.84 


142 ± 25=" 


NGC 4051 


0.00234 


0.191 ± 0.078 


0.0085 ± 0.0010 


87 ± 5 


81 


11.13 ± 0.25 




33.69 ± 0.52 


235 ± 32 


NGC 4151 


0.00332 


4.57 ± 0.52 


0.076 ± 0.040 


95 + 7 


181 


77.64 ± 1.55 


139 ± 26'" 


243.5 ± 1.5 


169 ± 13''' 


NGC 4593 


0.00900 


0.98 ± 0.21 


0.123 ± 0.040 


129 ± 15 


109 


3.48 ± 0.35 




13.79 ± 0.86 




NGC 5548 


0.01717 


6.54 ± 0.26 


0.324 ± 0.324 


192 ± 15 


174 


3.23 ± 0.21 


200 ± 33 


12.80 ± 0.66 




NGC 7469 


0.01632 


1.22 ± 0.14 


0.525 ± 0.024 


142 ± 3 


153 


12.07 ± 0.63 




40.63 ± 3.06 


179 ± 22* 


PG 0003+199 


0.02578 


1.42 ± 0.37 


0.724 ± 0.067 




119 


1.21 ± 0.30 


253 ± 50 


6.44 ± 0.28 




PG 0026-1-129 


0.142 


39.3 ± 9.6 


10.47 ± 1.45 






0.38 ± 0.12 




2.71 ± 0.38 


291 ± 51 


PG 0804-1-761 


0.1 


69.3 ± 8.3 


8.710 ± 1.604 






< 0.40 




1.93 ± 0.25 


308 ± 57 


PG 0844-1-349 


0.064 


9.24 ± 3.81 


2.239 ± 0.206 






< 0.72 


267 ± 56 


1.47 ± 0.21 




PG 0953-1-414 


0.2341 


27.6 ± 5.9 


16.60 ± 2.29 




231 


< 0.48 




< 1.14 




PG 1211-(-143 


0.0809 


14.6 ± 4.4 


5.623 ± 0.777 






0.53 ± 0.14 


252 ± 58 


1.45 ± 0.26 


209 ± 40 


PG 1226-1-023 


0.1583 


88.6 ± 18.7 


91.20 ± 10.50 




321 


2.61 ± 0.42 


504 ± 64 


9.17 ± 0.33 


423 ± 61 


PG 1229-1-204 


0.0603 


7.32 ± 3.52 


1.202 ± 0.138 


162 ± 32 




0.93 ± 0.07 




2.59 ± 0.45 




PG 1307-1-085 


0.155 


44.0 ± 12.3 


7.586 ± 0.700 




219 


0.67 ± 0.06 


262 ± 38 


1.96 ± 0.33 




PG 1411+442 


0.0896 


44.3 ± 14.6 


4.266 ± 0.393 




175 


0.64 ± 0.09 


284 ± 54 


1.87 ± 0.24 


228 ± 52 


PG 1426+015 


0.0865 


129.8 ± 38.5 


5.248 ± 0.846 


185 ± 67 


197 


1.20 ± 0.17 


292 ± 46 


4.37 ± 0.56 


250 ± 38 


PG 1613+658 


0.129 


27.9 ± 12.9 


9.550 ± 1.099 




230 


1.16 ± 0.21 


202 ± 43 


5.28 ± 0.25 


206 ± 34 


PG 1617+175 


0.11244 


59.4 ± 13.8 


3.020 ± 0.556 


183 ± 47 




< 0.77 




< 1.62 




PG 1700+518 


0.292 


78.1 ± 17.4 


42.66 ± 2.95 






< 0.39 




< 1.82 




PG 2130+099 


0.063 


45.7 ± 5.5 


2.884 ± 0.266 


172 ± 46 




4.06 ± 0.20 


318 ± 44 


10.00 ± 0.31 





Note. — The M^jj and -Lopt values are from Peterson et al. (2004), Bentz et al. (2006; 2007), and Denney et al. (2006). The ct* values were extracted from the Call 
triplet (Onken et al. 2004; Nelson et al. 2004) in the Seyferts and from the CO band heads (Dasyra et al. 2007) in the QSOs. The ctjq jjjj values are from Peterson et 
al. (2004), Shields et al. (2003; 2006), and Bonning et al. (2005). For sources with two measurements of ct, the average value was used. The MIR line detection limit was 
set to 3 times the noise standard deviation. 

^The velocity dispersions of these lines were extracted from SWS spectra presented in Sturm et al. (2002). 
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